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Abstract 

I address the question of the fluctuations in fishery landings. Using the fishery 
statistics time-series collected by the Food and Agriculture Organization of the 
United Nations since the early 1950s, I here analyze fishing activities and find two 
scaling features of capture fisheries production: (i) the standard deviation of growth 
rate of the domestically landed catches decays as a power-law function of country 
landings with an exponent of value 0.15; (ii) the average number of fishers in a 
country scales to the 0.7 power of country landings. I show how these socio-ecological 
patterns may be related, yielding a scaling relation between these exponents. The 
predicted scaling relation implies that the width of the annual per capita growth- 
rate distribution scales to the 0.2 power of country landings, i.e. annual fluctuations 
in per capita landed catches increase with increased per capita catches in highly 
producing countries. Beside the scaling behavior, I report that fluctuations in the 
annual domestic landings have increased in the last 30 years, while the mean of the 
annual growth rate declined significantly after 1972. 

Key words: scaling law, growth dynamics, capture fisheries, stock exploitation, 
Fishstat Plus, FAO 



1 Introduction 



Fluctuations in fishery la ndings have been a collective human concern ((Botsford et al 



19971 IPaiily et~aTH 2002b Variation in supply (stock size) is a central problem 



in determining the optimal organization of the fishing industry and in how gov- 
ernments can best regulate fishing effort. Since the industrialization of fishing, 
fisheries scientists have been subject to intense political pressures: "Determine 
the causes of fluctuations in fish catches." For example, in the mid-1800s the 
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Norwegian parliament asked zoologists Michael Sars and George Ossian Sars 
to investigate the biology of Norwegian fisheries and to answer why the cod 
catches from the Lofoten Islands in Northern Norway fluctuated so greatly 
f)Smithl . ll99l . 



Nonlinear dynamics have been applied to understanding the ups and downs of 
populations. In the mid-1920s, an Italian biologist Umberto D'Ancona found 
the increase in predator fish and decrease in prey fish of various species in the 
Adriatic Sea during the World War I period from a statistical study of the 
numbers of e ach species sold on the fi sh markets of three ports, Fiume, Trieste, 
and Venice ( D'Ancona . 19261 1954 ). A famous Italian mathematician Vito 
Volterra, asked by D'Ancona whether he could come up with a mathematical 
model of predator-prey population dy namics, develop ed the first and simplest 
model of predator-prey interactions Q. Because an American 

mathematical biologist Alfred J. Lotka produced t he two-specie s interaction 
model independently and at about the same time (|LotkaLll925h . it is called 
Lotka- Volterra model. The nonlinear model predicts a cyclical relationship 
between predator and prey numbers. 



Thus, the jagged oscillation in populations was nothing new to ecologist, 
but before the 1970s, they put most of the patterns down to the stochas- 
tic, unaccountable effects: year-to year fluctuations in fish stocks are usu- 
ally attributed to variability in recruitment, competition, predation, ch anges 
i n fishing e ffort, and other sources of so-called environmental noise. iMavl 
f 1974 Il976l ) discovered that a density-dependent (i.e. nonlinear), determinis- 



tic model of ecology could pr oduce complex, random-looking patterns "chaos" . 
Conklin and Kolberd (J1994J) researched the effects of market demand on hal- 
ibut populations of Pacific Ocean fisheri es, and provided a mod el of chaotic 



dynamics for the Pacific halibut fishery. IWeisbuch et al.l (|l997l ) studied the 



influence of biological and socio-economic mechanisms on the amplitude of 
resource depletion, and showed not only that the depletion is most often due 
to inertial factors in the dynamics of exploitation, i.e. capital and labor over- 
shoot, but also that capital and labor inertia results into dangerous oscillations 
of the resource level, eventually finishing with a fishery crash. 



A central question in ecology for last decades has been how to quantify the 
relative impo rtance of stochastic and nonlinear factors for fluctuatio ns in pop- 
ulation size (jSaether et all 1200(1 iTurchinl . 120031: iHsieh et all bOOfih . Charac- 
terizing the fluctuation patterns is the major challenge for time series analysis 
in fisheries and for modeling the process of fish production. One major ap- 
proach to dealing with ecological complexity is to reduce the system to one 
or a few species: a number of detailed analysis of specific systems have cap- 
tured the patterns of fluctuat ions in fish stocks and catches, including noise- 
driven fluctuations and chaos (IWilson et al. L 1991 L Matsuda et al L 199ll 1992 ; 



Takeuchi et Z[ll99i logins et alHl997UBi0rnstad et allll99g[ 



Dixon et al 



1999: iBiOrnstad and Qrenfeli l200lh . 



Inspired by work of Keitt and Stanley ( 1998| ) and Keitt et al. ( 2002f ) , I here ex- 
plore an alternative approach to socio-ecological system dynamics, i.e. popula- 
tion dynamics within which fisheries are embedded. From the macro-ecological 
viewpoint, I analyze the time-series of annual fishery landings (all commercial 
species items fished inclusive) of 244 countries in total during the period 1950- 
2002, after excluding non-fish vertebrates such as whales. The data on annual 
domestically- landed catch are taken from Fishstat Plus ((FAQ! . l2004al ) . a set 
of statistical database compiled by the Food and Agriculture Organization 
(FAO) of the United Nations. FAO is the only one institution that maintains 
global fisheries statistics. The country-structured time-series of landed catch 
is a complex source of information. Annual landed catch is the aggregate of 
a large number of complex interactions between fishers and the stock. Some 
aspects of landed catch are biological phenomena (biomass of fish vulnerable 
to the gear) and technical efficiencies (use of machinery power and modern 
equipments such as echo sounders and fish detection devices); others are eco- 
nomic (market price of fish and the cost of applying effort), legal (regulations 
applied to fishing effort) and some are behavioral (skill of the fisher targeting 
the stock and accurate completion of catch records). Besides the fluctuations 
in fishery landings, I investigate the employment statistic s based on the FAO 
data on the number of commercial and subsistence fishers ( FAO . 19991 2004b| ). 
I report two power-law relationships between the following quantities for the 
capture fisheries sector: 



(i) the variability of landings time-series and the country landings 

(ii) the number of fishers in a country and the country landings 



In order to understand the scaling structure for the stock exploitation, a rela- 
tion between th ese power-law ex ponents is hypothesized, such as the Widom 
scaling relation ( Widom . 1965a fj). 



In the paper, symbols x and n are used for the annual domestically-landed 
catch and the number of fishers in countries, respectively. 



2 World fisheries trend 



Fisheries have recently become a topic fo r the media with global audiences , 
as well as f or the scientifi c comm unity [e.g. iBeddington and Kirkwood (2005); 
Ehrhardtl (|2005h : ISvmesI (2005)], with a focus on the wretched history, the 



poor current practice, and the future scena rios. For instanc e , the New Orleans 
Times-Picayune's eight-part 1996 series bv lMcQuaid et al.1 (J1996J) on "Oceans 
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of Trouble" won the Pulitzer Prize for public service in 1997, in reporting on 
threats to the world's fisheries. 



In 1950, the newly founded Food and Agriculture Organization of the United 
Nations began collection of global statistics. Fisheries in the early 1950s were 
at the onset of a period of extremely rapid growth. Everywhere that industrial- 
scale fishing (mainly trawling, but also purse seining and long-lining) was in- 
troduced, it competed with small-scale, or artisanal fisheries. Throughout the 
1950s and 1960s, this huge increase of global fishing effort led to an increase in 
catches so rapid that their trend exceeded human population growth, encour- 
aging an entire generation of fisheries managers and politicians t o believe that 



launc hing more boats would automatically lead to higher catches f)Paulv et al. 
20021 ). 



In 2002, fisheries provided direct employment to about 38 million people and 
accounte d for more th an 15 percent of the total human consumption of animal 
protein (|FAOL l2004bf) . Globally, first-sale capture-fishery revenues produced 
about 78 billion US dollars. Total world trade of fish products increased to an 
export value of 58.2 billion US dollars, continuing the last decade's underlying 
4.5 percent annual growth. The world's population has been increasing more 
quickly than the total food fish supply; as a result the average per capita 
fish supply has declined. Total production from marine and inland capture 
fisheries for the world reached 93.2 million tonnes (84.5 million tonnes marine 
and 8.7 million tonnes inland) [excluding China, 76.6 million tonnes from 
capture fisheries (70.1 million tonnes marine and 6.5 million tonnes inland)]. 



Figure 1A shows the trend of the global harvests with capture fisheries, demon- 
strating that the global capture production has been exponentially growing. 
The annual growth rate of global catches (defined as R = ln[X(t + 1)/X(t)}, 
where X(t) and X(t + 1) are the total capture production in the years t 
and t + 1) changed discontinuously at the critical time 1972, correspondent 
to the collapse of the Peruvian anchoveta fishery, signaling an abrupt transi- 
tion of bio-economic regime. The fisheries exploitation has expanded in two- 
exponential mode, and possibly reached a plateau with no transient period 
contrary to a logistic growth. 



In line with the increase in fisheries production, over the last three decades 
employment in fisheries and aquaculture has continued to increase in many 
countries. The world number of fishers and fish farmers has been exponentially 
(i.e. in a Malthusian mode) growing at annual rate of 3.6 percent (excluding 
China, 2.8 percent) for the period 1970-2002, depicted in Fig. IB. The world 
fishers population has been growing at a faster rate than the earth's human 
population (see inserted figure of Fig. IB). Since 1970 to 2000, the ratio of 
fishers and fish farmers to the whole population has increased from 0.332 to 
0.588 percent (excluding China, 0.349 to 0.490 percent). 
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Since there is the possibility of China's misreporting ()Watson and Paulvl .12001). 
for further analysis, I use the data with the exclusion of China's reported 
catches (including China does not change the results). 



3 Scaling of country growth rate 



In Fig. 2, shown are the probability distributions of the annual domestic land- 
ings (A) and the number of people actively fishing and fish farming in coun- 
tries (B), both which are consistent with a log-normal distribution 



P{x) 



X\ 2-ksI 



exp 



ln 2 (x/x) 
2sl 



where x and s x denote a geometric mean and a geometric standard devia- 
tion of the distribution, respectively. The log-normal relationship remained 
unchanged for the time period studied. This log-normal distributio n of the 
count ry landings is suggestive of Gibrat's law on market structure ((Gibratl 
193lh . It is interesting to note that, contrary to city sizes (resident population) 
or US firm sizes (employees and revenue), the country population engaged in 
fishing and fish farming does not reduce to the Zipf 's power-law distr i butio n 
characterized by the absence of a characteristic size |Zrol Il949t fAxtelll . l200lh . 
However, the conditional probability distributions of fishers population for 
countries with the same landings show "fat tails" with characteristic scales 
(for details, consult §4). 

In the limit of small annual changes in x, we can define the relative growth 
rate as 



= In 



x(t + 1) 
, x(t) 



x(t+l) -x(t) 
xjt) 



(2) 



where x(t + l) and x(t) are the annual landings of a country in two consecutive 
years. I calculate the conditional distribution p(r\x) of annual growth rate for 
countries with landings x. The distribution p(r\x) of the growth rate from 
1950 to 1971 is shown in Fig.3A, depicting a simple tent-shaped form on a 
semi-logarithmic scale. The distribution is not Gaussian, as expected from the 
Gibrat model, but rather is exponential: 



p(r\x) 



y/2, 



a r \x) 



exp 



V2\r- 



a r \x 




(3) 



where r(x) and a r (x) are the average growth rate and the standard deviation 
of r for countries with landings equal to x, respectively. An implication of this 
result is that the distribution of the growth rate has much broader tails than 
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would be expected for a Gaussian distribution. Notice that the fluctuation ay 
in growth rate r is independent of x. 

In Fig.3B, shown is the empirical conditional probability density of r for coun- 
tries from 1972 to 2002 with approximately the same landings, suggesting that 
p(r\x) has the same functional form, consistent with the Laplace distribution 
[Eq.(3)], with different widths <J r (x), for all x. In these years, contrary to the 
period of 1950-1971, we expect that the statistical properties of the growth 
rate r depend on landed catch x: the growth-rate distribution depends on the 
initial value of x. 

The study of the fluctuations has been shown to give important information 
regarding the underlying processes responsible for the observed macroscopic 
behavior. I therefore analyze the fluctuations in the growth rate of the country 
landings. While different plots for the values of landings fall onto one another 
for the 1950-1971 period, there is clear dependence of oy on the value of landed 
catch for both the periods of 1972-1989 and 1990-2002. Figure 4A shows 
that for the periods 1972-2002, a r (x) statistically decreases with increasing 
x. The width of the growth-rate distribution, cr r (x), scales over three order of 
magnitude, from landings about 10 3 tonnes up to landings about 10 6 tonnes, 
as a power law: 

a r (x) = Ax'? (4) 

where A 0.27 and (3 & 0.15 for the both investigated periods (landed value 
of catch unit in thousand tonnes). 

What is remarkable about Eqs.(3) and (4) is that they govern the growth rates 
of a diverse set of countries. Fishery labor productivity and capital intensity 
vary widely among countries. They range not only in their amount of landings 
but also in what and where they catch: the harvest included in the database is 
classified according to more than 1000 commercial species items, and according 
to the inland or marine area. Indeed, the data of growth rate for a wide range 
of production values plotted in terms of scaled coordinates collapse onto a 
single curve (Fig.4B), if the ordinate and the abscissa are chosen as 

p sc = V2a r (x)p(r\x) and r sc = V2 [r - f(x)] /cr r {x), (5) 

respectively. Accordingly, the probability density function of growth rate r for 
a country with x tonnes landings, detrended by the average growth rate f(x), 
should scale as 

MrV)d/ = /r (^)d^, (6) 

where the scaling function f r is exponential, and the detrended growth rate is 
denoted by r' [= r — r(x)]. 
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3.1 At-year growth rate 



We now question whether Eq.(3) of landed catches in countries is still valid 
for the growth rates for longer time scales. Consider the At-year growth rate 



TAt 



:ln 



In 



x(t + At) 

x(t+iy 



(7) 



x(t) 



+ In 



x(t + 2) 
x(t + l) 



+ • • • + In 



x(t + At) 
x(t + At- 1) 



Is then ta* for large At Gaussian? Figure 5A shows the distributions of growth 
rates for time scales from 1 year up to 30 years over the entire period (1950- 
2002) for all countries (excluding China). In Fig.5B shown are the probability 
distributions of the normalized growth rate defined as 



TAt ~ TAt 

9 At = , 

0~At 



(8) 



where TaI and a At are the average and the standard deviation of the At-year 
growth rate. We see that the distributions of gAt are, to a good approximation, 
symmetric with respect to zero. Figure 5C shows the positive tails of the com- 
plementary cumulative distributions P(> gAt) °f the normalized growth rates 
for the same time scales, compared to the standardized Gaussian distribution 
with a mean of and a standard deviation of 1. Because of the central limit 
theorem, the convergence to a Gaussian distribution on longer time scales is 
observed, as the Gibrat model predicts. 

We can further test the convergence to Gaussian behavior by analyzing the 
moments of the distribution of normalized growth rates. A more global quan- 
titative measure of distributions is in terms of the moments: 



= [\9At 



(9) 



where (. . .) denotes the average. The fractional moments of the normalized 
growth rates for At = 1, 10, 20 and 30 years (Fig.5D) show clear indication of 
convergence to Gaussian behavior with increasing At. However, it is remark- 
able that they do not satisfyingly approach the Gaussian moments even for 
large At = 30 years. 

Moreover, notice that the distribution of annual (At = 1) growth rates r\ for 
all countries and all years can be well-described by an exponential (Fig.5A), 



p[r\) 



y/2, 



o-i 



cxp 



V2\r- 



ri 



(10) 
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with rT = 0.04 and o\ = 0.27, while the conditional probability density of 
the yearly growth rates for countries with approximately the same landings x 
depe nds on x [Eqs.(3) and (4), Figs. 3 and 4]. Using a saddle point approxima- 
tion ()Lee et al. . 1998f ). we may integrate the conditional Laplace distribution, 
Eq.(3), over a log-normal distribution of the country landings, Eq.(l), and 
obtain Eq.(10). 



4 Scaling structure of the fisheries sector 



Let us analyze employment data in the fisheries sector in order to gain fur- 
ther insight into the organization of fisheries structure of countries. Figure 6A, 
showing a scatter plot of the number of fishers versus country landings, sug- 
gests that the conditional probability density p{n\x) to find a country with 
total catch x landed by n fishers has a long right tail. I here propose a model 
for the conditi onal distributions of fishers p opulation in terms of "stretched 
exponentials" ( Laherrere and Sornettel . 19981 ): 







" / n\ a ' 


) a 


— exp 




n 


Vn / 


\n J . 



'ID 



such shat the complementary cumulative distribution is 
P(> n\x) = I p(n'\x)dn f 



exp 



n 
n 



(12) 



where n denotes a reference scale of the distribution, expected to depend 
on the landings, no = no(x). Stretched exponentials are characterized by an 
exponent a smaller than one. The meaning "reference" is that from no all 
moments can be determined, e.g. the means of n and n 2 are given by 



n 



no a 



(13) 



and 



n- 



2n^- 1 r(2a" 1 ), (14) 

respectively, where T(q) is the gamma function [equal to (q — 1)! for q inte- 
ger]. When a is small, n will be much larger than n . The expression (12) 
of the stretched exponential complementary CDF means that the stretched 
exponential distribution is qualified by a straight line when plotting the neg- 
ative log of complementary CDF, — lnP(> function of n on a 
double-logarithmic scale. As shown in Fig.6B, lines appear: the values of the 
exponent a range from 0.60 to 0.86, numerically evaluated from Eqs.(13) and 
(14). Accordingly, the fishers population for coun tries with the s ame landings 
follows a stretched exponential distribution [aka IWeibullI (|l95lh distribution 
with exponent lass than one]. 
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Moreover Fig.6A suggests that the typical (i.e. average) number of fishers 
increases as a power law with the country landings x. We make the hypothesis 
that the conditional probability density p(n\x) obey the scaling relation 

P(»|*)dn = U (£) d (£) , (15) 

where f n denotes a scaling function. To test the scaling hypothesis in Eq.(15), 
we plot the scaled quantities 

p sc = x a p[n\x) vs n sc = n/x a , (16) 

and observe the data collapse with a exponent a = 0.71 ± 0.01 (Fig. 7). Since 
the data exhibit scaling, we identify a universal scaling function fitted with a 
stretched exponential distribution 

fn(n 9c ) = an*- 1 exp (-/£) . (17) 

Notice that the average number of fishers scales with country landings x as 

n(x) = / np{n\x)dn cx x a . (18) 
J o 

As shown in Fig.6A this is indeed the case. 



The multiplicative cascade model ( Frisch and SornetteL 1997 ) may explain the 



stretched exponential distribution of fishers: the probability distribution func- 
tion of the product of m independent identically distributed positive random 
variables with an exponential- like distribution leads to stretched exponentials 
[see also Laherrere and Sornettd ( 19981 )]. The reciprocal number of levels in 



the multiplicative cascade reads the index a = mT x . Interpreted within the 
multiplicative cascade model, the value of the exponent a ~ 0.7 corresponds 
to about one to three generation levels in the generation of a typical fisher. 



Notice that the conditional distribution p{n\x) reduces to a stretched exponen- 
tial distribution, while the distribution of the country population engaged in 
fishing and fish farming is log-normal (Fig.2B). It is numerically demonstrated 
(Fig. 8) that a log-normal distribution of the number of fishers approximately 
results from a convolution of a stretched exponential distribution p{n\x) and 
a log-normal distribution of the country landings, Eq.(l): 



Pin) 



p{n\x)P(x)dx 



oo n - 

a — 

o nf, 



a-l 



exp 



n 
n 



X\ lixsi 



■ exp 



ln 2 (x/x) 



2*2 



dx (19) 



with n oc x a . 
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5 Modeling scaling behavior in the exploitation dynamics 



Total catch of a country is made up of individual fisher activity in the country. 
Fisheries management ranges along a continuum from community to govern- 
ment as sole managers of the resource, and the field of impact of their decisions 
can be local, regional, provincial, or national. Fishers get involved in their re- 
spective management arrangements in a country. 



We now consider the following two limiting cases. 



(i) There are strong correlations between the fishers in a country, following that 
the growth dynamics are indistinguishable from the dynamics of structure- 
less organizations. As a result, there is no country-size dependence of the 
growth fluctuation a r , following (3 — 0. 

(ii) The annual landed catches by individual fishers change independently of 
one another. Then, if the number of fishers in a country were proportional 
to its domestic fishery product, the growth fluctuations as a function of 
domestic ally-landed catch wou ld decay as a power law with an exponent 



(3 = 0.5 ((Buldvrev et all Il997l ). This is because the standard deviation of 



the sum of n independent quantities grows like ^fri (i.e. catches by n fishers). 



As to the stock status during the period 1 950-1971, there was much room 
for continued expansion of capture fisheries ( Grainger and GarciaL 199fi| ). The 
fishing policy those days would have been virtually that "Strive for as much 
fish as possible from the resources." All the decisions in a country would have 
been perfectly coordinated as if they had been all dictated by a single "boss" . 
The state of stock exploitation is then expected to have corresponded to the 
case (i): annual catches by individual fishers in a country were correlated 
strongly each other. 



In 1972 or subsequent years, the observed decrease in fluctuations with coun- 
try landings, Eq.(4), is considerably slower (i.e. (3 ~ 0.15) than in the case (ii). 
Our empirical results suggest an important consequence for dynamics of stock 
exploitation: although countries with large fisheries tend to diversify into a 
wider range of fishing activities leading to smaller relative fluctuations in coun- 
try landings, fisheries catches reported are more variable than what would be 
expected if the number of fishers would increase linearly with the country 
landings — which would correspond to (3 — 0.5. 
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5.1 Scaling relation between exponents 



The country landings does not scale in a simple linear fashion with increas- 
ing number of fishers, but instead the typical number of fishers scales with 
landings as Eq.(18). The total capture production is broken down by num- 
bers of landing sites (i.e. fishing ports) in a country. The fishery products in 
landing sites are landed by numbers of local fishers. These production equals 
the marketing volume of fish at local landing sites. Let us analyze fishery 
statistics of Japan on the marketing of fishery products in major landing are as 



((Ministry of Agriculture Forestry and Fisheries of Japanl . 12002. 2003. 2004 ). It 



is virtually apparent in Fig. 9, showing the probability distribution of catches 
landed at 203 ports (A) and the probability density of annual growth rate 
of port landings (B), that the same behavior holds for subsystems, i.e. the 
single country's fisheries sector, as the world capture fisheries sector. The sub- 
systems or hierarchical systems of capture fisheries sector are expected to be 
characterized by the same fundamental laws. 

To account for the increased variability for highly producing countries (i.e. 
(3 ~ 0.15 for years 1972-2002), we suppose that in a country with x tonnes 
domestic fisheries production, the width of the annual per capita growth-rate 
distribution, denoted by o~o(x), scales with landings x as 

o-o (x) oc x 1 (20) 

with an exponent 7. Motivated by scaling behavior in the exploitation dynam- 
ics , I expect that t here exist scaling relations among three exponents a, (3 and 
7 (|Widoml . Il965a 3). Assume that the per capita growth rate is independent 



(uncorrelated) , identically distributed, though individual fishers may not be 
right fishing units that produce catches independently of one another; e.g. 
there are 132 thous and fishery establishments while 238 thousand fishers ex 



ist in Japan in 2003 f Mi nistry of Agriculture Forestry and Fisheries of Japan . 



20051 ). The central limit theorem is then applied to give 



a r [x) = A_ , (21; 
nix) 



leading to a scaling relation 



0=?-7- (22) 



For a = 0.71 ±0.01 and (3 = 0.15 ±0.02, Eq.(22) predicts 7 = 0.21 ±0.03. It is 
remarkable that, contrary to country landings (Fig.4A), the per capita landings 
of highly producing countries fluctuate more than do those of low producing 
countries [Eq.(20) with 7 > 0]: fishers in a country with larger landings risk 
more than do fishers in a country with smaller landings in capture fisheries. 
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6 Discussion 



Studying the biological properties of the FAO global fisheries statistics on 
landed c atches has given ecological insight to the st a te of world fish stocks: e.g . 
in Refs flPaulv et all Il99fit ICaddv and GaribaldvL |200(A IPaulv et all l2005h . 



global conclusions drawn state that a general decline in mean trophic level of 
marine and inland landings has occurred in many regions, reflecting a change in 
catch composition from long-lived, high trophic level, piscivorous bottom fishes 
toward short-lived, low trophic level invertebrates and planktivorous pelagic 
fishes. In this paper, from the socio-ecological viewpoint, I have analyzed the 
FAO data and found the universal scaling behavior in the dynamics of stock 
exploitation. Moreover, we empirically find from Fig. 6 and Eq.(18), that the 
per capita production of countries with landings x scales as 

oc X 1 "", (23) 



n[x) 



showing the increased per capita catches for highly producing countries (i.e. 
1 — a > 0). However, fishing business may have a high-risk premium: the 
larger their expected revenues in highly producing countries, the larger the 
fluctuations in annual per capita revenues, which is predicted by Eq.(20) with 
7 > 0; whereas the empirical data demonstrate decreasing fluctuations in 
annual domestically-landed catches with increased country landings [Fig.4A 
and Eq.(4) with (3 > 0]. The expected fishing risk (i.e. fluctuations in the 
individual catches) in a country with x tonnes annual domestic landings scales 
with the individual expected catches: a (x) oc [x/n(x)} s , with an exponent 
5 = 0.7 ± 0.1, which seems to agree with the value of the exponent a. 

Though the growth era of the capture fisheries production has been ended 
during the last decade, the world fishers have been growing in number (35.8 
percent from 1990 to 2002). The fishers population to total world popula- 
tion ratio has been increasing for the past thirty years and stood at 60.7 
fishers per one thousand people in 2002. In most developing countries of low 
and middle-income, the number of people employed in the fisheries sector has 
been growing steadily. In industrialized economies offering occupational al- 
ternatives, the numbers of fishers have been on a declining trend or at best 
stationary. For instance in two important fishing countries, Japan and Norway, 
fishers have decreased in number by 34.3 and 19.7 percent, respectively, during 
the 1990 to 2002 period. The increasing trend of world fishers population may 
suggest that wages received by fishers (mostly in developing countries) exceed 
the level of wages that fishers could expect to earn in alternative employment 
opportunities, as commonly understood in economic analysis from the static 
viewpoint of the common-property fishery. Because fishers would be making 
money more than their opportunity costs (i.e. resource rent is positive), addi- 
tional fishers would be attracted to the fishery, until it reaches the "bionomic" 
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equilibrium (jGordonl Il954h where total sustained revenue equals total sus- 
tained cost. Contrary, if net revenues from fishing are negative, fishers leave 
the fishery. 

The exploitation dynamics of fish stocks exhibits two-phase behavior, which 
may be according to stock status together with management policy. There was 
a change in scaling behavior of the standard deviation of production growth 
rates at 1972. It may be more than coincidence that the emergence of the 
United Nations' Convention on the Law of the Sea, in the late 1970s, which 
enabled countries to claim exclusive economic zones reaching 200 nautical 
miles into the open sea, put the responsibility for fisheries resource manage- 
ment squarely with maritime co untries, thus endi ng many decades of fighting 



over traditional fishing grounds ((Johnstonel 11977). 



Another possible cause for the observed transition at 1972 in scaling struc- 
ture is that, having depleted large valuable stocks, fishing has redirected 
some effort and added a lot of it on other species lower down the food web. 
The strategy was advocated in the 1970s to increase fisheries production 
f Sprague and Arnold . 19721 Garcia and Grainger . 2005f ). The ph enomenon of 



chang e in catch composition, i.e. "fishing down marine food webs" flPaulv et al. 
1998), c ontrasts str i kingly with D'Aancona's observation on Mediterranean 



fisheries ([D'AnconaL I1926T ) that the cessation of the fishery during World War 
I from 1914-1918 changed the species composition in favor of fish-eating fishes 
(the frequency of predators decreased with an increase in fishing after the war). 
The decline in mean t rophic level of the world fish stocks leads to increasing 
variability in catches ( Paulv et all I2OO2T 1 . Indeed, we find from the analysis 
of the FAO data for all countries (excluding China) that the country landings 
fluctuate more for the periods 1972-2002 than do those for the periods 1950- 
1971. The standard deviation of growth rate of the annual domestically-landed 
catches for 30 one-year periods between 1972 and 2002 shows an increase of 
more than 10 percent over that for 21 one- year periods between 1950 and 
1971: 

OY,i95o-i97i = 0.138 ± 0.004, cf ril972 _ 2002 = 0.156 ± 0.003, (24) 
whereas the average yearly growth rate declined significantly after 1972: 



0.077, 



0.026. 



(25) 



Besides highly variable recruitment of planktivorous pelagic fishes that is 
strongly influenced by environmental fluctuations, the schooling behavior com- 
mon among such pelagics is a possible cause of variation in catch in purse 
seine fisheries. Some fishers catch more than others: the distribution of fish- 
ers' yearly total lan ded value of ca tch in a commercial purse seine fishery is 
skewed to the right ( Hilborn . 1985J ). School-size distributions of pelagic fishes, 
i.e. catch-per-haul distributions in purse seine fisheries, are generally heavy- 
tailed, and a typical school-size (rigorously defined as a crossover size from 
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power-law to exponential decay of the school-size dis tribution) is proportional 
to the abundance of the fish stock ( Niwa . 20031 2004 ). seeming to be consistent 
with Eq.(20) with 7 > 0. 

The global potential harvest (e.g. maximum sustainable yield) may be one 
of th e causes for the observed tr ansition in scaling behavior of fishery land- 
ings. Grainger and Garcia (|l996| ) estimated the most likely potential of con- 
ventional marine species to range from 80 to 100 million tonnes, which has 
indeed been reac hed probably in the 1970s and is unlikely to change in the 
next 20-30 years ( Garcia and Grainger . 2005| ). Undeveloped resource fisheries, 
producing much less than their potential, dec reased rapidly to zero by the mid- 
dle of the 19 70s (iG rainger anHGarciai lSh. The latest analysis of the state 
of resources ( Garcia et al. . 2005| ) indicates that, in 2003, approximately half 
the world's stocks are exploited at or close to their maximum, and about 25 
percent of them are exploited either below or above such maximum. 



The same scal ing laws as found h ere describe a number of other complex 
organ izations ( Stanley et al. . 20021). such as b usiness firms ( Stanley et al. 



_996), gross domest ic product (|Lee et all Il998l). university r esearch budgets 
( Plerou et al. . 1999j ). and bird populations (|Keitt et al. . 2002| ). suggesting the 
underlying universal mechanism. It is no wonder that commercial fisheries 
world-wide resemble competitive economic activities. In the same way that 
market forces operate in business, the stock exploitation appears to keep com- 
petition among fishers strong enough to leave ca pture fisheri es with a shrinking 
resource base. According to FAO press release ( FAQ . 1998J ). most small-scale 



fishers find it increasingly difficult to survive in an over-exploited environment. 



Amaral et al. ( 1998J ) proposed a model for the growth of organizations with 
complex internal structure. In the model, each organization (e.g. firm, uni- 
versity, or bird population) is made up of units. The units composing the 
system have a complex evolving structure, e.g. business firms competing in an 
economy are comprised of divisions; universities are composed of schools or 
colleges; population of breeding birds living in a given ecosystem is composed 
of groups living in different areas. Their model postulates that the units grow 
through an independent, Gaussian-distributed, random m ultiplicative process. 



Later this postulate empirical ly proved not to be ac curate (M atia et all 12004). 
Thus the model proposed by Amaral et al. (1998) can now be regarded as a 
first step towards the explanation of the growth dynamics of organizations. 
Scaling concepts offer an avenue to study heterogeneous assemblies in ecolog- 
ical and economic systems for which the microscopic proces ses are (probably) 
not k nowable, except in terms of their statistical properties ([Chave and Levin . 
20031 ). 
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Figure legends 



Fig.l World fisheries trend including (o) /excluding (•) China. (A) Semi-logarithmic 
scale plot of global capture trend. Exponentially growing behavior with two 
growth-rates is observed. There is a knee in the curve at 1972, probably fig- 
uring the boundary between two distinct regimes of the bio-economies. The 
solid lines are fits of an exponential function [oc exp(i?£)] to the data (•) 
with average yearly growth rates R = 0.059 ± 0.001 and 0.024 ± 0.001 
for the periods 1950-1971 and 1972-1989, respectively. A plateau in the 
1990s follows immediately (to be precise R = —0.002 ± 0.003). Fitting the 
data (o) gives R = 0.058 ± 0.001 and 0.017 ± 0.001 for the periods 1950- 
1971 and 1972-2002, respectively (depicted by broken lines). (B) Trend 
of the world fisher population N. The ordinate represents log 10 N. The 
broken and solid lines depict exponential fits to the data including and 
excluding China, giving population growth rate Rn = 0.035 ± 0.002 and 
Rn = 0.028 ± 0.001 for the period 1970-2002, respectively (annual growth 
rate of the fisher population world-wide is defined as same as the global 
harvest case). The inserted figure shows the fisher-to-population ratio. The 
earth's hu man population data from 1970 to 2000 we re retrieved from the 
web-site o f lUnited Nations Population Division! (|2004h . The human popula- 



tion has been growing at an average rate of 1.67 percent per year (excluding 
China, 1.74 percent) since 1970 to 2000. 

Fig.2 Fisheries production and employment statistics. Data exclude China's statis- 
tics. (A) Semi-logarithmic scale plot of probability distribution of country 
landings, P(x). The abscissa is chosen as logarithmically scaled landings 
ln(x/x)/s x . The data from 6 one- year periods between 1950 and 2000 are 
plotted, where geometric means x of landings (in thousand tonnes) and geo- 
metric standard deviations s x are as follows: x = 6.0, 9.8, 15.4, 16.9, 20.1 and 
20.4, and s x = 2.5, 2.6, 2.7, 2.9, 2.9 and 2.9, from years 1950 (□), 1960 (0), 
1970 (A), 1980 (■), 1990 (♦) and 2000 (▲), respectively. Log-scaled land- 
ings fall onto the standardized Gaussian distribution (dashed line) with zero 
mean and unit variance. (B) Semi-logarithmic scale plot of probability dis- 
tributions of the number of fishers and fish farmers in countries, P(n). The 
abscissa is chosen as logarithmically scaled fishers population \n(n/n)/s n , 
where geometric means n of the number of fishers (in thousand fishers) 
and geometric standard deviations s n follows: h = 6.7, 8.1 and 10.2, 

and s n = 2.1, 2.1 and 2.0, from years 1970 (A), 1980 (■) and 1990 (♦), 
respectively. Log-scaled fishers population fall onto the standardized Gaus- 
sian distribution (dashed line). The log- normal distribution holds robustly 
for fishers as well as landings. 

Fig. 3 Exploitation dynamics of fish stocks. (A) Semi-logarithmic scale plot of 
probability density p(r\x) of the annual growth rate for eight different bins 
of fisheries production (in tonnes): x < 10 2 ' 5 (o), 10 25 < x < 10 3 (□), 



19 



10 3 < x < 10 3 - 5 (0), 10 3 - 5 < x < 10 4 (A), 10 4 < x < 10 4 - 5 (■), 10 4 - 5 < 
x < 10 5 (♦), 10 5 < x < 10 5 - 5 (A), 10 5 - 5 < x (•). The data are aggregated 
over all 21 one-year periods between 1950 and 1971. The distribution decays 
with fatter tails than for a Gaussian. The broken line is a fit of a Laplace 
distribution [Eq.(3)] for the period, with r = 0.077 and a r = 0.138 ± 0.004. 
It is depicted that p(r\x) does not depend on x. (B) Semi-logarithmic scale 
plot of probability density of the growth rate for countries with small and 
large landings in the years 1972-1989/1990-2002: 10 3 < x < 10 4 (0/ A) 
and 10 5 < x < 10 6 (♦/A) in tonnes. The broken and solid lines are Lapla- 
cian fits using the means r = 0.035 and 0.016 and standard deviations 
o r = 0.203 ± 0.004 and 0.117 ± 0.004, calculated from data for countries 
with small and large landings in the years 1972-1989, respectively. Fits with 
Eq.(3) for small and large landings in 1990-2002 give the means r = 0.034 
and 0.005 and standard deviations a r = 0.219 ± 0.009 and 0.118 ± 0.007, 
overlapping broken and solid lines in the graph. 

Fig. 4 Scaling in stock exploitation dynamics. (A) Standard deviation of the one- 
year growth rates of the 21 annual intervals from the period 1950-1971 (o), 
of the 17 annual intervals from the period 1972-1989 (•), and of the 12 
annual intervals from the period 1990-2002 (■), as a function of the landings 
x (plotted on a double-logarithmic scale). Each data point indicates the 
width a r (x) of the growth-rate distribution for countries with approximately 
the same landings (binned evenly spaced on a logarithmic scale). The long- 
dashed horizontal line depicts the average of the growth fluctuations for 
years 1950-1971 (cr r = 0.138), showing f3 = 0. The solid and broken lines 
are fits with Eq.(4) to the data from years 1972-1989 {(5 = 0.14 ±0.02) and 
1990-2002 (P = 0.15 ± 0.02), respectively. There are cutoffs at the small 
landings in power laws. Error bars of one standard deviation about each data 
point are shown. (B) Scaled probability density p sc as a function of the scaled 
growth rate r sc for seven different bins of production in the years 1972-1989 
(in tonnes): 10 3 < x < 10 35 (□), 10 35 < x < 10 4 (0), 10 4 < x < 10 4 - 5 (A), 
10 4 - 5 < x < 10 5 (■), 10 5 < x < 10 5 - 5 (♦), 10 5 - 5 < x < 10 6 (A), 10 6 < x (•). 
The values are scaled by using the measured values of r(x) and a r (x), and 
plotted on a semi- logarithmic scale. All the data collapse upon the universal 
curve p sc = exp (— |r sc |) (broken line). A data collapse can be obtained for 
1990-2002 as well. 

Fig. 5 Growth rates of country landings for longer time scales. The data points 
are aggregated over the entire period, 1950-2002, for all countries (exclud- 
ing China). (A) Probability densities p(r& t ) of At-year growth rates are 
shown for At = 1 (□), 10 ((>), 20 (a) and 30 years (■), plotted on a semi- 
logarithmic scale. Averages and standard deviations of the At-year growth 
rate are as follows: rXi = 0.04, 0.39, 0.73 and 1.08, and a At = 0.27, 0.80, 0.99 
and 1.11 for the same time scales. The solid line indicates a Laplace distribu- 
tion [Eq.(10)] with rj~ and a\. It is visually apparent that the distributions 
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approach a Gaussian as At increases. (B) Semi-logarithmic scale plot of 
probability distributions of the normalized growth rates g^t- The broken 
line shows the standardized Gaussian distribution. (C) Double-logarithmic 
scale plot of complementary cumulative distributions P(> g& t ) for the pos- 
itive tails of the probability densities shown in Fig.B. The complementary 
CDFs for the negative tails have the same functional forms as for the pos- 
itive tails. (D) Fractional moments for the normalized growth rates for the 
same time scales: At = 1 (long-dashed line), 10 (medium-dashed line), 20 
(short-dashed line) and 30 years (dotted line). The solid line shows the 
Gaussian moments. 

Fig. 6 Statistical analysis of the number of fishers. Data exclude China's statistics. 

(A) Number of fishers versus country landings (in tonnes). The ordinate and 
the abscissa are chosen as log 10 n and log 10 x, respectively. Data are aggre- 
gated for seven years (1990-1996). The rectangle in gray reads the inter- 
quantile range of binned data and the slit indicates the median, showing that 
the conditional distribution p{n\x) is skewed to the right, where the bins are 
chosen equally spaced on a logarithmic scale as x G [10 fe ~ a5 , 10 fc+0 - 5 ) with 
k — 3, . . . , 6. The dashed line shows a power-law fit to scatter plot data, sug- 
gesting that the typical number of fishers scales with country landings as x a . 
From the slope of the dashed line, the exponent a is evaluated at 0.68±0.08. 

(B) Stretched exponential law in the conditional probability of n given x, 
p(n\x). The negative logs of complementary CDF, — lnP(> n\x), are plot- 
ted versus decimal log of the number of fishers, log 10 n, on a semi-logarithmic 
scale. The data points fall onto straight lines. The countries are partitioned 
into four groups according to landings (in tonnes): 10 25 < x < 10 35 (()), 
10 3 - 5 < x < 10 4 - 5 (A), 10 4 - 5 < x < 10 5 - 5 (♦), 10 5 - 5 < x < 10 6 - 5 (A). The 
solid lines are stretched exponential fits to the data, representing power-law 
curves, -lnP(> n\x) = (n/n ) a , where a = 0.86 (0), 0.72 (A), 0.66 (♦) 
and 0.60 (A), estimated from Eqs.(13) and (14). The estimated reference 
scales no are depicted in Fig. A with symbol (o), which closely follow the 
dashed line. The exponent characterizing stretched exponential distribu- 
tions should be the exponent a in the power-low dependence of the average 
number n(x) of fishers in a country with annual domestic landings x as 
suggested in Fig. A (dashed line). A power-law fit with Uq(x) = Bx a to es- 
timated n 's (solid line in Fig.A) gives a = 0.71 ± 0.01 and B = 11 ± 2, 
exhibiting a good agreement with values of the exponent estimated from 
Eqs.(13) and (14). 

Fig.7 Scaling in employment statistics for the period 1990-1996. (A) To test the 
scaling hypothesis in Eq.(15), the scaled probability distributions p sc (= 
Bx a p) are plotted versus scaled number n sc (= n/Bx a ) of fishers (double- 
logarithmic scale). The data for four different bins are scaled by using the 
values of landings x = 10 3 (0), 10 4 (a), 10 5 (♦), 10 6 tonnes (A), with 
a = 0.71 and 5 = 11. The conditional PDF p{n\x) exhibits scaling, i.e. the 
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data collapse onto a single curve, yielding a scaling function f n . The solid 
line depicts the stretched exponential distribution [Eq.(17)]. (B) Shown is 
the data collapse of the four curves in Fig.6B. The negative logs of comple- 
mentary CDF, — lnP(> n sc ), are plotted versus scaled number n sc of fish- 
ers (double-logarithmic scale). The solid line depicts the power-law curve, 
-lnP(>n sc )=< c . 

Fig. 8 Recovering a log- normal distribution from a stretched exponential. Semi- 
logarithmic scale plot of probability distribution of the number of fishers 
for all countries (excluding China). The abscissa is chosen as log 10 n. The 
data points are the average over the period 1990-1996, and the dashed line 
is a fit of the log-normal PDF to the data with geometric mean n = 12.9 (in 
thousand fishers) and geometric standard deviation s n = 2.09. Numerically 
integrating the conditional PDF p(n\x) using stretched exponential function 
Eq.(17) with a = 0.71 over a PDF of the country landings using Eq.(l), 
a log-normal distribution of the number of fishers is approximately recov- 
ered (solid line). The reference size of the stretched exponential distribution 
p(n\x) is n (x) = Bx a with B — 11. The distribution of the landing data 
for years 1990-1996 has the geometric mean x — 17.1 (in thousand tonnes) 
and the geometric standard deviation s x = 3.22. 

Fig. 9 Statistics on the marketing of fishery products in landing sites (203 ports) 
in Japan. (A) Semi-logarithmic scale plot of probability distribution of local 
capture production £. The abscissa is chosen as logarithmically scaled land- 
ings In (£/£)/«£, where geometric means £ of landings (in thousand tonnes) 
and geometric standard deviations are as follows: £ = 4.01, 3.60 and 3.33, 
and s e = 1.75, 1.70 and 1.66, from years 2000 (□), 2001 (■) and 2002 (<)), 
respectively. Log-scaled landings fall onto the standardized Gaussian distri- 
bution (dashed line). The port landing data exhibit the log-normal prob- 
ability distribution, same as the distribution of country landings (Fig. 2 A). 
(B) Semi-logarithmic scale plot of probability density function of annual 
growth rate rg of port landings. The points are the average over the years 
2000-2002. We find the same tent shape distribution as the distribution of 
country landings growth rate plotted on a semi-logarithmic scale. The solid 
line is an exponential fit, p(r^) = (y^og) -1 exp V^\r^ — r^/a^J, with av- 
erage growth rate = —0.06 and standard deviation = 0.30 ±0.02. The 
dashed line is a Gaussian fit with standard deviation <r Gauss = 0.49. 
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Figures 

Fig.l World fisheries trend including (o) /excluding (•) China. 

Fig. 2 Fisheries production and employment statistics. 

Fig. 3 Exploitation dynamics of fish stocks. 

Fig.4 Scaling in stock exploitation dynamics. 

Fig. 5 Growth rates of country landings for longer time scales. 

Fig. 6 Statistical analysis of the number of fishers. 

Fig.7 Scaling in employment statistics for the period 1990-1996. 

Fig. 8 Recovering a log-normal distribution from a stretched exponential. 

Fig. 9 Statistics on the marketing of fishery products in landing sites in Japan. 
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